Switching exponent scaling near bifurcation points for non-Gaussian noise 
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We study noise-induced switching of a system close to bifurcation parameter values where the 
number of stable states changes. For non-Gaussian noise, the switching exponent, which gives 
the logarithm of the switching rate, displays a non-power-law dependence on the distance to the 
bifurcation point. This dependence is found for Poisson noise. Even weak additional Gaussian noise 
dominates switching sufficiently close to the bifurcation point, leading to a crossover in the behavior 
of the switching exponent. 
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Physical systems display generic features near bifurca- 
tion parameter values where the number of stable states 
changes. In this range the dynamics is controlled by a 
slow variable, a soft mode. Its noise-induced fluctua- 
tions are comparatively large. They ultimately lead to 
switching of the system from the stable state. Close to 
a bifurcation point the switching rate W becomes appre- 
ciable even where far from this point it is exceedingly 
small, for a given noise level (for example, for given tem- 
perature). The high sensitivity of the rate to the system 
parameters has been broadly used to determine parame- 
ters of Josephson junctions and Josephson junction based 
systems Jll-Qj nanomagnets [B43i mechanical nanores- 
onators [J, and recently in quantum measurements 

The analysis of switching conventionally relies on the 
assumption that the underlying noise is Gaussian. Then 
the switching exponent Q, i.e., the exponent in the ex- 
pression for the switching rate W ex. exp(— Q), displays a 
power-law dependence on the distance to the bifurcation 
point in the parameter space t], Q (x rj^ [l|, [l^, For 
systems in thermal equilibrium ^ = 3/2 for a saddle- node 
bifurcation and ^ = 2 for a pitchfork bifurcation. This 
applies also to systems far from equilibrium [13 . 

Recently, there has been much interest in large fluctu- 
ations and switching induced by no n- Gaussian noise 
[T9| . Such switching can be used to determine the noise 
statistics [20l - l26j . However, the features of the switching 
rate near bifurcation points have not been explored. Yet, 
one may expect that the 77-dependence of the switching 
exponent will differ from that for a Gaussian noise and 
will be very sensitive to the noise statistics. 

In this paper we study the behavior of the switching 
exponent Q for systems driven by Poisson noise. Such 
noise is often encountered in photon statistics and in the 
statistics of current in mesoscopic conductors. We show 



that the scaling is not described by a simple power law, 
and the overall 77 dependence of Q is much weaker than 
for Gaussian noise. Surprisingly, if in addition to Poisson 
noise the system is driven even by a comparatively weak 
Gaussian noise, sufficiently close to the bifurcation point 
this noise dominates and there occurs a crossover to the 
standard scaling of Q for Gaussian noise. 

Generally, one would expect that, unless it is very 
weak, a Poisson noise would make a stronger effect on 
the switching rate than a Gaussian noise. This is so, 
since switching is a rare event on the scale of the char- 
acteristic relaxation time of the system tr, it requires a 
large fluctuation, whose probability is determined by the 
tail of the noise distribution. Such a tail is less steep for 
a Poisson noise than for a Gaussian noise. 

The "takeover" of the switching rate by a weak Gaus- 
sian noise close to a bifurcation point is a more subtle 
effect. It emerges because of the qualitatively different 
ways the fluctuations leading to switching occur for Gaus- 
sian and Poisson noises. This can be understood from the 
equation of motion for the slow variable q, 



q^-U'{q) + f{t). 



(1) 



Here, U{q) is the effective potential; for the saddle- 
node and pitchfork bifurcations, U = U^'^'^\q) and U = 
U^P^^q), respectively [13], with 
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m' + ^q^. (2) 
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For ?7 > the system has a stable state qa at the mini- 
mum of U (q) (or two such states, for a pitchfork bifurca- 
tion) and a saddle point qs at the local maximum of U (g); 
for 77 = these states merge together. The results can be 
immediately extended also to the case [/^p^^ —U^p^\ 
the subcritical pitchfork bifurcation where for ?; = a 
stable state merges with two unstable states. 

The force f{t) in Eq. ([T]) is noise. We will consider 
the cases where f{t) is a Poisson noise, f{t) = fp{t) 
with fp(t) = g'^n^i^ ~ ^")' ^ white Gaussian noise. 
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fit) = fait), {fG{t)fG{t')) = 2DS{t - i'), or a combina- 
tion f{t) = fp{t) + fcit)- Both g and D are assumed 
small, so that the switching rate W ^ — U"{qa)- 
The noises are (5-correlated in time, because the system 
motion is slow; they are also independent of q, since of 
interest is a small region in the system phase space [3] ■ 
We start with the qualitative picture of Poisson-noise 
induced switching. Here switching occurs only for the ap- 
propriate pulse polarity, (qs ~ qa)/g > 0- A single noise 
pulse shifts the coordinate q by pulse area g. Switching 
requires driving q from to q^, from where the sys- 
tem will switch with probability ^ 1/2. The necessary 
'^o ~ (<75 — 9a)/5 pulses should occur within time < tr, 
so that the system cannot relax back to the attractor be- 
tween the pulses. The probability of such pulse sequence 
is equal to (vtr)"" cyi'p{~i>tr)/nol, where v is the average 
pulse frequency. By construction, this is the probabil- 
ity to switch in time tr-, it is ~ Wtr, and thus gives 
the switching rate. The corresponding estimate of the 
switching exponent Q = Qp for no 3> 1, 1'tr is 

Qp ~ [{qs - qa)/9] In [{qs - qa)/9Vtr\ . (3) 

From Eq. \qs — qa[ k 7?^^^, whereas tr ^ i]^^^ with 
^ = 3/2 and ^ = 2 for the saddle- node and pitchfork 
bifurcations. It is seen from Eq. ^ that Qp oc rj^^'^ but 
contains a large ?7-dependent logarithmic factor; the ratio 
{is ~ qa)/gi^tr is the large parameter of the theory. 

White Gaussian noise fG{t), on the other hand, leads 
to switching by providing a force that overcomes the de- 
terministic force U'{q). From Eq. ©, [U'{q)\ < rf--^!"^ 
for qa < q < qs- Since the probability of noise realization 
IS oc exp 

[-/ dtf^{t)/AD] ^ and the duration of the 
needed noise outburst is ^ tr, by setting this probability 
to be ^ Wtr and fa ^ U' , we obtain for the switching 
exponent Q ~ Qc the familiar expression [ij, 

Qg = Cv^/D (4) 

(C — 4/3 and C = 1/4 for the saddle-node and pitchfork 
bifurcation, respectively, see below). 

With decreasing 77, the distance between the stationary 
states qs — qa oc rj^^^ decreases slower than the determin- 
istic force, \U'{q)\ < 77^^^/^. Therefore an outburst of 
Poisson noise required for a transition decreases slower 
than that of Gaussian noise. Respectively, as seen from 
Eqs. Q and Q, the switching exponent decreases much 
slower for Poisson noise than for Gaussian noise. As a 
result, for small Gaussian-noise intensity D ^ g, the 
switching rate is determined by Gaussian noise for suffi- 
ciently small 7], whereas for larger -q it is determined by 
Poisson noise. 

The crossover between Gaussian- and Poisson-noise 
dominated switching is clearly seen in Fig. [1] It presents 
the results on switching near a saddle-node bifurcation 
point in the presence of both noises. Close to the bi- 
furcation point the slope of In Q vs. In 77 is 3/2, as for 
Gaussian noise. However, as rj increases the slope ap- 
proaches that for purely Poisson-noise. 
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FIG. 1: Color online. The crossover of the switching exponent 
Q from the Poisson-noise to Gaussian-noise dominated behav- 
ior with the decreasing distance 77 to the saddle-node bifurca- 
tion point. The solid line shows the result of Eqs. © and ((T)) 
for the Gaussian noise intensity D = 0.009, the average Pois- 
son pulse frequency 1/ = 0.09, and the Poisson pulse amplitude 
g = 0.185; the data points are the results of numerical simula- 
tions in which Q was found from W using the Gaussian-noise 
prefactor, W ~ {r}^^^ /n) exp{—Q). The dashed and dotted 
lines were obtained from Eqs. ^ and ((7]) by setting D = or 
g = 0, respectively. 

We now give a quantitative theory of the switching ex- 
ponent and find its 77-depcndence for Poisson noise as 
well as study the crossover between the Poisson- and 
Gaussian-noise asymptotic. As a starting point, we use 
a generalized Fokker-Planck equation for the probability 
density of the system, which follows from Eq. ([T]) 

dtp = d, [U'{q)p{q)] + Dd'-p{q) + v [p{q - g) - p{q)] . (5) 

The last two terms in this equation describe the effect 
of the Poisson noise, i.e., of uncorrelated pulses with av- 
erage frequency v that shift the coordinate by g. These 
terms have the same form as reaction terms in the mas- 
ter equation for a reaction system, with q and g being 
the number of species and the change of this number in 
a reaction, respectively [29j . except that in the present 
case q is continuous; the analogy applies where g is small 
compared to the typical scale of q, in particular compared 
to qs — qa- In this case a Poisson noise in the Langevin 
equation mimics reactions in reaction systems. 

For small D and g, the switching rate is determined 
by the probability current away from the initially occu- 
pied attraction basin [s^l- The current is independent 
of time for ■C t -C and, to logarithmic accu- 

racy, is given by the quasistationary probability distri- 
bution at the saddle point. To find this distribution, 
we solve Eq. ([5]) in the eikonal approximation: we set 
p{q) = exp[— s(g)], assume that s{q) — s{qa) ^ 1, and 
keep the leading order terms in s [for example, we disre- 
gard dgS compared to {dgs)'^]. This leads to the equation 
for s{q) of the form H{q, dqs) = 0, where 

H{q, p) = ~pU'{q) + Dp'-v[l~ exp(5p)] . (6) 
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Equation ([6|) maps the problem of the quasistationary 
distribution of the fluctuating system onto the problem 
of Hamiltonian dynamics of an auxiliary conservative sys- 
tem with coordinate q, momentum p, and Hamiltonian 
H{q,p). The switching exponent is Q = s{qs) ~ s{qa), or 



p{q)dq. 



(7) 



Here, qa,Qs a-re the stationary states shifted to allow 
for the nonzero mean of the Poisson noise, qi — qi « 
vg/U'^qi) with i = a,S. The momentum p{q) in Eq. ([7]) 
is the nontrivial solution of equation H{q,p) = 0, it gives 
the trajectory of the auxiliary system which goes from 
9a to qs. Equations (O could be obtained also by 
finding the probability density of the most probable re- 
alization of noise fjt) necessary to drive the system from 
to qs 0, Hll , albeit such calculation would be 
somewhat more involved. 

For purely Poisson noise, i.e., where D — 0, for 
small pulse area g, equation H(q,p) = gives p « 
g^^{lnv{q) + \n[\nv{q)]}, where v{q) = U'{q)/vg. This 
estimate applies provided v{q) 3> 1; we disregarded 
higher-order corrections oc l/ln[z;(g)]. The condition 
v[q) ^ 1 holds in much of the region between and 
qs except the immediate vicinities of qa,QSj i-G-, for 
\Q — Qa\, \q — qs\ ^ >^\g\/U"{qa)- This follows from the es- 
timate |L'''(g)| < ri^~^/^ in the central part of the interval 
iqa,qs) and the inequality {qs-qa)/g > i^U = v/U"{qa) 
discussed above. 

Keeping in p{q) the leading order term and replac- 
ing U'{q) by its maximal value between qa and qs, from 
Eq. ^ we obtain the following estimates for the switch- 
ing exponent for the saddle-node and pitchfork bifurca- 
tion, respectively. 



(8) 



The parameters k^''"-', k'p^-' in the arguments of the loga- 
rithms are ^ 1. A simple choice k^™^ = 2 and k^p*^ = 1 
gives a close agreement of Eq. ^ with the results ob- 
tained by numerically solving equation H ~ and 
Eq. ([71), which are shown in Figs. [2] and [3] For the cho- 
sen ly, the difference is < 10% for Q > 20 and 5 > 0.1 . 
Equations ^ justify the estimate Eq. ([3]). 

The power-law factor in Eq s. (Et is the same for 
both types of the bifurcation points j32| . It is determined 
simply by the distance between qs and qa- However, the 
arguments of the logarithms are different. The logarith- 
mic factors significantly change the switching exponent 
compared to a simple power-law scaling Qp oc rj^^^/g. 
For example, for the pitchfork bifurcation for g = 0.2 in 
Fig. |3] the logarithmic factor varies from 3.6 for r/ = 0.5 
to 2.2 for T] = 0.2. 

The results of the asymptotic theory were compared 
with numerical simulations of the system Eq. ([T]). As seen 



12r 



-|(sn) 




0.1 



0.2 



FIG. 2: Color online. The switching exponent Qp"' for a 
Poisson noise as a function of the distance rj to the bifurcation 
point for the saddle-node bifurcation. The mean frequency of 
noise pulses is f = 0.1. The solid lines are obtained from 
Eq. ((7|) using a numerical solution of equation H{q,p) — 0. 
The data points show the results of numerical simulations of 
switching; the plotted quantity is ln{ri^^^ /Wn). 



12 



-i(pf) 







g=0.05 


* 


g=o.i 


□ 


g=0.15 


A 


g=o.2 





g=0.25 



0.25 



0.5 



0.75 



FIG. 3: Color online. The switching exponent Qp*' for a 
Poisson noise as a function of the distance rj to the bifurca- 
tion point for the pitchfork bifurcation. The mean frequency 
of noise pulses is = 0.05. The solid lines are obtained 
from a numerical solution of equation H{q,p) = followed 
by integration, Eq. ([7]). The data points show the results of 
numerical simulations of switching; the plotted quantity is 
\n{ri/WTvV2). 



from Figs. [2] and m they are in excellent agreement. In 
determining the switching exponent from the switching 
rate VF, we used the pref actor in W which coincides with 
that for white Gaussian noise [1^. We note that, for 
the supercritical pitchfork bifurcation described by the 
potential J7'-P^-'(g) in Eq. a unipolar (all pulses of the 
same sign) Poisson noise leads to switching from only 
one of the two coexisting stable states at the minima of 
U^^^\q). For the subcritical bifurcation described by the 
potential —U'^^^^q), Poisson noise always leads to decay 
of a metastable state. 

In the opposite case of purely Gaussian noise, g — 0, 
the nontrivial solution of equation H(q,p) = is p = 
U'{q)/D. It leads to a power-law dependence of the 
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switching exponent on rj described by Eq. (|4]) and to the 
appropriate values of the constant C in this equation. 

Equation ([6]) describes quantitatively the crossover 
from Poisson- to Gaussian-noise dominated switching as 
T] approaches the bifurcation value rj = 0. If the Gaussian 
noise is much weaker than the Poisson noise, D <^ g, far 
from qa,qs the momentum p(q) is close to the Poisson- 
noise (D ~ 0) solution provided rj is not too small, so that 
M{r]) > D/g, with M(?/) = / In [7]'^-'^^'^ /i^g) (we 

use that, near its maximum, \U'\ ^ 77^^^/^). On the other 
hand, where M{r]) <ti D/g the momentum is determined 
by the Gaussian-noise {g = 0) solution. The position of 
the crossover on the 77-axis is given by M{r]) ^ D/g. Not 
surprisingly, for such 77 the switching exponents for the 
purely Poisson and purely Gaussian noises become of the 
same order of magnitude, Qp ~ Qg- This argument is 
confirmed by the data in Fig. [TJ 

The situation where a weak Gaussian noise is present 
even where other sources of noise are dominating is typ- 
ical for practically any physical system. Such noise very 
often comes simply from the coupling of the system to 
a reservoir that leads to energy dissipation. The results 
of the paper explain why near bifurcation points there is 



often observed the power-law scaling typical for Gaussian 
noise even where this noise is comparatively weak. 

In conclusion, we have considered noise-induced 
switching due to a non-Gaussian noise near two generic 
types of the bifurcation points: saddle-node and pitch- 
fork. In contrast to the case of Gaussian noise, where 
the switching exponent scales as a power of the distance 
to the bifurcation point 77, for a non-Gaussian noise the 
exponent generally displays a more complicated depen- 
dence on 77. We have found it for a Poisson noise, in 
which case, along with a power-law factor, the exponent 
has a large logarithmic factor. It turned out that even a 
weak additional Gaussian noise becomes the major cause 
of switching sufRciently close to the bifurcation point. A 
qualitative and quantitative description of the crossover 
from Poisson to Gaussian noise controlled switching and 
of the r^-dependence of the switching exponent are in full 
agreement with numerical simulations. 
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